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Abstract 

We study limits for the mass varying neutrino model, using constrains from super- 
nova neutrinos placed by the r-process condition, < 0.5. Also, we use this model in 



' a supernova environment to study the regions of survival probability in the oscillation 

. space parameter (tan^ 9 and Attiq), considering the channel — t- and Ue — i^s- We 

observe that r-process nucleosynthesis does not put severe constraints in our model, 
although new regions of isocurves of probability appear when our parameterization 
. includes an inversion in the hierarchy between the neutrino eigenvalues rrii and m^. 

o 

m 

^ '. 1 Introduction 
(N 

^ ' Evidences from experimental data of type la supernovae (SNIa) [H [2], cosmic microwave 

background (CMB) radiation [3l H] and large scale structure (LSS) [5] point out that our 
universe is in accelerated expansion. One possible explanation for this dark energy issue 
^ ! (70% of the energy content of our universe) would be a cosmological constant (A), or simply, 
^ I the non-zero vacuum energy, which is quite close to the critical cosmological energy density, 
Pc ~ Pa ~ 4 X 10~^'^GeV^. However one critical point of A is the necessity of a fine tuning, 
once that according to theoretical expectations, p^ac is 10^° — 10^^° larger than the magnitude 
allowed in cosmology. In reference [6], Dolgov emphasizes: 

1. Why vacuum energy, which must stay constant in the course of cosmological evolution, 
or dark energy, which should evolve with time quite differently from the normal matter, 
have similar magnitude just today, all being close to the value of the critical energy 
density? 
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2. // universe acceleration is induced by something which is different from vacuum energy, 
then what kind of field or object creates the observed cosmological behavior? Or could 
it be a modification of gravitational interactions at cosmologically large distances ? 

From the questions above, other candidates for dark energy have been proposed, such as: 
quintessence, K-essence, tachyon field, phantom field, dilatonic dark energy, Chaplygin gas 
(for excellent reviews see [Tj.jll. 

Neutrinos are powerful tools for astronomy investigation ^ and they play a crucial role 
in the mechanism of supernova explosion, according to the standard scenario [9]. Also, 
due to their peculiar features, supernova environments are potentially interesting to obtain 
information about possible new physics. 

Based on the fact that the scale of the neutrino mass squared difference ((0.01 eV)^) 
is similar to the dark energy scale, the possibility of neutrinos coupling with a scalar field 
naturally arises. In this picture, the total energy of the fiuid can vary slowly as the neutrino 
density decreases [10]. Consequently, the coupling with a scalar field causes a variation 
on the neutrino mass (dynamical field) with the neutrino density and perhaps with the 
baryonic matter [H]. This picture has been called Mass Varying Neutrino models, MaVaN. 
Bjalde et al. [12] showed that the proposed MaVaN models, in the non-relativistic regime, 
can be threatened by the strong growth of hydrodynamic perturbations associated with a 
negative adiabatic sound speed squared, many articles have treated neutrino oscillation via 
the mechanism of MaVaN [l3l[TIl[T5l[T6l[l71[T8l[T9l[20l|2T]. We observe, nevertheless, 
that these articles focus on solar and atmospheric neutrinos, as well as observations from 
accelerators and reactors. However, the assumption of MaVaN mechanism in the neutrino 
evolution in supernovae environments has not yet been fully explored. In fact, the only 
article that connects MaVaNs and supernova is Li et al. [22] which uses data from a type la 
supernova to limit the interactions between neutrinos and scalar dark energy. 

In this paper, therefore, we investigate the consequences of MaVaN mechanisms for the 
neutrino propagation in a supernova environment. We take into consideration that recent 
results from WMAP-7 [23] indicate that the number of species of relativistic neutrinos is equal 
to 4.34±0.87 which, in contrast with the standard model with 3 active neutrinos, suggests that 
there may be another species of neutrinos, for example, a sterile one Ug, possibly associated 
with SU(2) singlets. These sterile neutrinos can present masses ranging from a few eV to 
values larger than the electroweak scale [21]. Therefore, a natural and potentially interesting 
investigation can arise from the analyzes of neutrino oscillations Ug — )■ and — )■ z/^ in 
supernovae under the assumption that MaVaN models are implemented. We evaluate the 
consequences of such neutrino oscillations on the signal of Ue and in terrestrial detectors, 
as well as on supernova heavy-element nucleosynthesis. 

The crucial point of our analysis is the fact that MaVaN model can directly affect the 
relevant parameter for neutrino oscillations, namely, the squared mass difference Am^. We 
propose a phenomenological MaVaN model that is more convenient to test using neutrinos 
from supernova. Such parameterization allows to fit the amplitude of the Am^ variation, the 
position where the variation begins and also if the masses will increase or decrease along the 

^ There are alternative scenarios to explain the acceleration, such as modified gravity, where occurs an 
introduction of quantum effects - higher curvature corrections - to the Einstein-Hilbert action. 
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neutrino trajectory. We use this parameterization to evaluate the constraints on the oscilla- 
tion space parameter sin^ 26 x Attiq (where 6 and Amg are the mixing angle and the squared 
mass difference in vacuum, respectively) coming from the r-process nucleosynthesis condition, 
Ye < 0.5, and by the limit for the average survival probability coming from SN1987A data, 
(Pee) < 0.5. Problems and criticisms to these MaVaN models were pointed by Peccei in 

We conclude that, with our MaVaN parameterization, significant modifications happen 
in the probability space parameter for anti-neutrinos when the mass modification is slow 
and for an inversion of hierarchy during the propagation. Also, for neutrinos, new regions 
of exclusion are also found and the r-process nucleosynthesis is not affected by the MaVaN 
assumption. 

This article is organized as follows: section |2] presents the neutrino oscillation mechanism 
and the assumptions and approximations used in our work. Section Elpresents our proposition 
of a new parameterization of MaVaN model. Section H] contains our results and related 
discussion. Finally section [5] summarizes our conclusions. 



2 Neutrino Oscillation 

Some comments about the assumptions for the neutrino oscillations we are interested in this 
article are in order. From the neutrino phenomenology we know the existence of two small 
scales, Am^j^ and two non-small mixing angles, 621 e 6*23 and one small angle 6'i3. We are 
interested in oscillations between active electronic into sterile neutrinos as well as oscillations 
between active electronic into sterile anti-neutrinos in a supernova environment. For this we 
will include a the new scale, A for neutrinos oscillations with sterile neutrinos. These 
new scale need to much bigger then the small scales for active neutrinos, Attiq = 1 — 
10^ eV^ > Am^j^ to be compatible with constrains for sterile neutrinos. For an analysis 
of sterile neutrino oscillations in cosmological, astrophysical and terrestrial media with various 
mixtures and a wider range of Am? see |26] . 

Such oscillations must be treated considering that both neutrino and anti-neutrino con- 
version can be resonantly enhanced in matter, the so-called Mikheyev-Smirnov-Wolfenstein 
(MSW) phenomena [27] . Nevertheless, for a convenient choice of the squared mass difference 
in vacuum, which we will denote Amg, we are allowed to neglect the conversion to sterile 
neutrinos in the inner core, for all values of the mixing angles we considered |28]. Further- 
more, this same choice of the Amp parameters allow neglecting v ~ v forward-scattering 
contribution to the weak potentials in first approximation after neutrinos have escaped from 
the inner core of the stain 

The equation which governs the flavor neutrino eigenstates evolution along their trajectory 
inside the supernova can be written as: 



. d 


' ^e(r) ' 








" ^e(r) ' 


dr 


^s(r) 




-0e(r) 







(1) 



where 

0e(r) = -^(±y(r)i?- Am2cos20), (2) 



^ Recently issues of self neutrino interactions, in the time we are considering (for the cooling stage), were 
discussed by A. Fricdland in [5^. 
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and 



^=^sin2^. (3) 



4:E 

Under the assumptions described above, the matter potential can be written as follows 



V{r) = 2V2GF 



(4) 



In Eq. ([2]), the signal + is for neutrino, whereas the — signal is for anti-neutrinos. Attlq is 
the squared mass difference between two neutrinos mass eigenstates in vacuum; Gp is the 
Fermi coupling constant; 6 is the vacuum mixing angle and iVe-(r), Ne+{r) and Nn{r) are, 
respectively, the number densities of electrons, positrons and neutrons. Note that neutrinos 
(anti-neutrinos) undergo in a resonance when 

V{r) = ±=^cos2e. (5) 
2E 

The matter potential can be rewritten as 
where is the electronic fraction 



iV,(r) + JV„(r) ' 

The spectrum of neutrino emission and it can be written as [30] : 

dN ^ L 

dE F{ri)T^ e'^/^-^ + l ^ ' 

In the expression above, Eq. ([8]), 77 is the pinching factor and Firf) = Jl^ dxx^ / {e^~^ + 1). 
For this spectrum {E)/T ^ 3.1514 + 0.1250?7 + 0.042977^ + 0(?7^). Typical values of t] are 
rji,^ ~ 2, rjp^ ~ 3 and 77^ ~ 1 [31]. The average energy of neutrino species is {E^J ^11 MeV, 
{EpJ ^ 16 MeV, {E^J ^ 25 MeV (x = /U,r). To calculate the probability in regions that 
are not adiabatic, we numericaly solve Eq. ([T]). And for regions that are adiabatic we use 
the analytical expression for the survival probability 0.5 + 0.5 cos 20 cos 2^0 (^0 is the matter 
mixing angle in the point of neutrino production). The calculation of average probability is 
done with a average over the spectrum of neutrino emission described in Eq. |8l 
We discuss our MaVaN phenomenological model in next session. 



3 Mass Varying Neutrino model 

We will adopt a phenomenological approach in modeling the MaVaN mechanismH. We 
parameterize the effects of new physics directly in Am^, since we are interested in neutri- 
nos oscillation and their implications on the oscillation probability and the r-process in the 

■^In |32] it was proposed also a phenomenological construction of a MaVaN model, but the analysis was 
made in the cosmological context. 
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supernova. We present a model in which Am\j^y^^ = Am? can increase or decrease as the 
neutrino propagates in the supernova medium, with the following parameterization: 

Am^ = Ami - 1 7 — (9) 

where rii, is the neutrino density in the supernova environment, and the parameters n[l and 
7] control the profile of neutrino mass dependence. 

The mass variation occurs around n^, and for << (n,^ » n^) the value of Am^ 
tends to the asymptotic value of Am^ {Am^ — 6). The parameter rj controls how fast the 
asymptotic values are achieved. 

The neutrino evolution with MAVAN is described by changing Am^ -^■ Am^ in Eqs. (ni2][5|^ 
for neutrinos and anti-neutrinos. 

Our parameterization is very versatile in reproducing the behavior of any specific models 
which includes MaVaN. For instance, the model presented in [15] is very well reproduced by 
choosing S = 10 eV^, t] = 1 and = 34500 cm~^. 

Clearly, if 5 = we have the "standard" scenario, in which neutrinos have a constant mass. 
According to the model, if 5 > Am^, then we can obtain Am^ < somewhere inside the 
supernova. This tell us that the hierarchy is being inverted in some point of the propagation. 

Figures [1] shows the evolution of Am^ with the supernova radius (r) by Eq. (Q for some 
parameters choices. 



4 Results and Discussion 

It is possible to constraint the oscillation parameters, such as the mixing angle and the 
mass-squared difference, from the SN1987A data [331 ISl |35]. One of them comes from the 
spectrum of the observed events, probably dominated by Ue- From the observed energy 
spectrum of SN1987A, the analysis from Kamiokande indicate that the temperature of i>e is 
lower than the expectation ^B]- The second constraint is based that the expanding envelope 
driven by thermal neutrino wind of the supernova is a possible site of heavy nuclei formation 
beyond iron (r-process nucleosynthesis). Since the Ue and Ue conversion into a sterile state 
(or even into an active state) are different, the fraction of neutrons to protons, determined by 
Ue + P ^ + n and + n ^ e~ + p, is modified by the oscillation mechanism in a radius of 
about a hundred kilometers. Related to this constraint, a very simple bound is Fg < 0.5, which 
can be extracted from Eq. ([7]) simply considering a very neutron- rich environment [HTJ [3H] • 
In [38], it is shown that Ye has a dependence on the neutrino oscillations probabilities by 



^ 1 + P,{E,J/P,{E,J ' ' 

given in Eq. (17) of |SH]), where Py and P^, are, respectively, the survival probabilities of anti- 
neutrinos and neutrinos. From this expression, we clearly see that it is possible to constraint 
the oscillation parameters from Ye < 0.5. Also in |36], it is also discussed another constraint 
based on the fact that the first or the second event of Kamiokande could be related with Ue- 
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Figure 1: The evolution of Am^ from Eq.l^ is shown for Attiq = 10 eV^ , 6 = ±10 eV^ and 
= 10^^, 10^^ and 34500 cm~^ for the solid, dashed and dotted lines, respectively. The values for 
the parameter r] are shown in the plot. 



We divide this section into two parts. The first one concerns the analysis of the oscillation 
of active electronic (anti)neutrinos to sterile (anti)neutrinos. In the second section we discuss 
the possible limits imposed by r-process nucleosynthesis in the oscillation space parameter. 
We emphasize that we are seeking changes caused by our MaVaN model, presented in section 

El 

4.1 (Anti)neutrinos 

Here we calculate the survival probability (Pee) for the Ue i^s channel of oscillation. For 
this purpose the calculation of this survival probability is done analytically for regions where 
the neutrino propagation is adiabatic. When there is a non-adiabatic region we calculate 
numerically using Eq. ([1]). We emphasize that the MaVaN model changes the regions where 
resonances occur, since there is a modification of Am^ along the neutrino path and these 
values enter in place of AttIq in Eq. ([5|). The electronic potential profile we have used is for 
a supernova 1 second after supernova bounce [tph > Is), where the size of the neutrinosphere 
is approximately equal to the core size and the supernova is in the stage where heavy nuclei 
start to form. Based on the data from SN1987A and neutrino oscillation, several articles, 
such as [38], discuss that the survival probability cannot be less than 0.5. In this way it is 
possible to find regions of exclusion in the oscillation space parameters (AmQ,tan^^). We 
explore here < 6 < n/A, called light side, and 7r/4 < 6 < 7r/2, called the dark side, which 
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can be interpreted for neutrino oscillation with < 6* < 7r/4. So our results can also be 
extended for the neutrino oscillation case, where dark side of anti-neutrinos will be light side 
for neutrinos and vice versa. 

We will find probabilities in the context of the MaVaN parameterization, given by Eq. (Q. 
We consider four sets of (77, n1) values for each 6 (in eV^) which will assume values of 2, -2, 
20, -20. They are: (0.5,1 x lO^^ cm-^), (0.8,1.5 x lO^o cm-^), (3.0,1 x 10^^ ^^-3) g^^^ 
(10,5 X 10^^ cm~^). All these sets ensure that Arh^ reaches the vacuum value inside the 
supernova. Also they are only illustrative, i. e., another values could be used for similar 
results. 

We start our analysis with the curve evolution set (77, n^) = (0.5, 1 x 10^^ cm^^) for 
6 = 20 eV^. In Fig. [2] we note the presence of a new area in the region of tan^ 6 ^ 10^^ — 10^ 
and around AttLq ^ 1 — 10 eV^. Above this region the result is close to the situation without 
MaVaN (5 = 0). In Figs. [2|4f5| the thicker curves represent the case without MaVaN. The 
presence of new regions of probability is explained by the fact that Am^ reaches negative 
values, i.e., rh^ < rhi, and this creates new regions of resonance, which can be seen in Fig. [3l 
where resonance points are represented by the crossings between neutrino potential curve and 
the resonances curves, which are calculated by Eq. We emphasize that, until r ^ 13.3 km 
(bottom part of Fig. [3]), for the situation without MaVaN and normal hierarchy [m^ > mi), 
the first region of the potential is an important region for resonances of anti-neutrinos for 
< 6 < it/A. The top part of Fig. [3] shows also the important region of resonances when 
Tr/A<e < tt/2. 

However, if Am^ is negative in some region inside the supernova, up to three resonances 
may occur for i^e- It is precisely the creation of new resonant points that will generate new 
regions in the parameter space for a given value of probability. 

We also note in Fig. |2] that for very similar values of Attiq and a single value of tan^ 6, 
there may be very different values of survival probability. Fig. E] shows the electronic potential 
{Ve, solid line) and the resonance curves for Amg = 1 eV^ (dashed line) and Amg ^ lA eV^ 
(dotted line), for a value of tan^^ ?a 2.5 x 10~^ and 6 = 2eV'^. These values are merely 
illustrative and other values could have been chosen. For Amg ^1.4 eV^ we notice that 
there are two consecutive resonances, keeping the survival probability practically constant. 
For Attt-o = leV^ the survival probability changes when the electronic anti-neutrino crosses 
the first resonance, as it is expected. After the electronic anti-neutrino propagates and 
crosses the second resonance, the probability fall is interrupted at the exact moment when 
the transition of hierarchy occurs from a negative value to the positive one, resulting in a 
different value of probability compared with the Amg ^ lA eV^ case. 

In Fig. m which represents the state of evolution (0.8, 1.5 x lO^'' cm~^) for S = 20 eV^, 
we see the reduction of the effect, since the change in hierarchy is brought further back and 
some regions of new resonances are lost. For 6 = 2 eV^, the effects are not so relevant, since 
the durability of the inversion of hierarchy is small. 

For sets of (77, n°) equal to (3.0, 1 x 10^^ cm~^) and (10, 5 x 10^^ cm~^) there is a great 
coincidence of the survival probabilities with the situation without MaVaN. This happens 
because the evolution of Am^ is so fast in both sets that it reaches rapidly Attiq and kills the 
possibility of creation of new resonances. Since the curves of MaVaN and 6 = practically 
match, we do not show them here. 

Now we present the results when 6 assumes a negative value. In Fig. |5] we show the 
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Figure 2: Isocurves of survival probability (Pee)- The dashed represents Pee = 0.1, the solid one 
Pee = 0.5 and the dotted curve Pee = 0.7. The curve evolution set (rj, n^) = (^O-S, 1 x 10^* cmT^) 
for (5 = 20 eV^. In this figure, the thicker curves represent the case without MaVaN (6 = 0). 



results for 6 = -20 eV^ and (r/, n°) = (0.5, 1 x 10^* cm since for these values we get 
the most significant modification. We do not show the other situations, because there is 
a complete match among the curves of MaVaN and without MaVaN. For 6 = — 20eV^ the 
effect is more significant, especially for values of small Am^, because the effect of a higher 6 
is more important in these regions, since the resonance position can be modified or simply 
removed by MaVaN effects. Here the change of hierarchy does not happen, then we do not 
expect any new resonances and, therefore, very different and unusual regions defined by the 
probability curves will disappear, since the resonance regions of anti-neutrinos are still for 
r < 13.3 km (light side) and r > 13.3 km (dark side), when S is negative. Again we see that 
as Am^ reaches faster the value of Amp, the probability curves will coincide for the MaVaN 
and the usual situation, 6 = 0. We do not show for 6 = —2 eV^ since the effect is very small 
and can be neglected. 

4.2 r-process 

The production of heavy nuclei in supernova is still an open problem. The region between 
the protoneutron star and the escaping shock wave a few seconds after the bounce may be 
a good site for this process with high entropy and an excess number of neutrons. From /3 
reactions, UeP ne~^ and UeU — pe~ , a modification in the n/p fraction could happen if one 
considers neutrino oscillation [371 EH] . We will analyze what happens with this condition in 
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Figure 3: The solid curve represents the electronic potential (Ve), calculated from Eq. Dotted 
and dashed curves are the resonance curves for Attiq ~ 1.4 eV^ and AmQ = 1 eV^, respectively, 
with 5 = 2eV'^ and (ry,n° j=(t). 5,1 x lO^*^ cm~^) for tan^ 6 « 2.5 x 10~^. When there is a cross 
between these curves with the potential one, we have a resonance point. 



the context of our MaVaN parameterization and for z/g — z/j, and i/g — >■ i^s- We limit our 
analysis of neutrino flavor conversion to about r < 50 km, i. e., a region where there is not 
any shock wave and has a higher "influence" on anti-neutrinos. This radius is approximately 
the region where the r-process nucleosynthesis happens - for a review in this subject see |10]. 

Since the analysis is limited to r ^ 50 km, when there is a slow evolution of Am^ to 
the vacuum value, i. e., for the sets (0.5,1 x 10^^ cm~^) and (0.8,1.5 x 10^° cm"^), no 
resonances will happen for electronic anti-neutrinos since there is an inversion of hierarchy 
considering 5 = 2, 20 eV^. However, for neutrinos, there are resonances, but all of them 
are adiabatic with different initial conditions, so, the r-process with the MAVAN mechanism 
always give Yf, < 0.5, independent of oscillation parameters. Even for 6 = —2, —20 eV^ 
where the initial conditions are modified, the Yg < 0.5 is also maintained. 

For the sets, (3.0, 1 x 10^^ cm"^) and (10, 5 x 10^^ cm"^), considering any 6 = 2, 20, 
2,-20 eV^, where Am^ reaches very fast Amg, we obtain approximately the same Fig. 7 
of reference [38], since, as we commented before, these situations are very similar to the 
situation without MaVaN (5 = 0). 

The evolution of Am^ represented by the dotted line in Fig. [T] and its implication in the 
r-process nucleosynthesis and oscillation probabilities were discussed in |41j . 
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Figure 4: Isocurves of survival probability (Pee)- The dashed represents Pee = 0.1, the solid one 
Pee = 0.5 and the dotted curve Pee = 0.7. The curve evolution set (i], n^) = (0.8, 1.5 x 10^° cm~^) 
for 6 = 20 eV^. In this figure, the thicker curves represent the case without MaVaN (6 = 0). 



5 Conclusions and Outlook 

In this paper we proposed a new phenomenological parameterization for the variation of 
the relevant neutrino oscillation parameters generally present in MaVaN models. This param- 
eterization has the advantage of being very versatile and easily adaptable to any theoretical 
and phenomenological approach derived from MaVan models. Considering the formulation 
of Eq. iQ, we have analyzed the modification in the neutrino survival probability for the 
channels Ue — ^ i^s and z/g — )■ Ug compared to the case 6 = 0, which corresponds to the absence 
of MaVaN. In general, we observe that new regions in the neutrino oscillation parameter 
space appear when 6 assumes positive values. In particular, we observe this behavior for 
5 = 20eV^, and (r/, n°) equals to (0.5, 1 x 10^^ cm-=^) and (0.8, 1.5 x 10^° cm"^). When 
the amplitude gets smaller, then the effects also will be smaller. On the contrary, when 6 
assumes negative values, no significant modification, if we compare with the positive case, 
happens. 

In both cases cited before, the r-process nucleosynthesis condition Ye < 0.5 is always 
guaranteed and, consequently, no significant new constraints to the oscillation parameter 
space appear from the enforcement of this condition. 

Furthermore, when the evolution is settled by (3.0, 1 x 10^^ cm^^) and (10, 5 x 10^^ 
cm^^) no significant modifications were introduced with respect to standard, where no Ma- 
VaN is present, since the evolution of Am^ quickly reaches Attiq, approximately mimicking 
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Figure 5: Isocurves of survival probability (Pee)- The dashed represents Pee = 0.1, the solid one 
Pee = 0.5 and the dotted curve Pee = 0.7. The curve evolution set (rj, n^) = (^O-S, 1 x 10^* cmT^) 
for 6 = —20 eV^. In this figure, the thicker curves represent the case without MaVaN (6 = 0). 



6 = probabilities. In the last section, \4.2\ we pointed out that these sets will reproduce the 
Fig. 7 of [3H] for the r-process. 

We point out that these kind of studies, without the MaVaN approach, were made in the 
past by |371l38]- The reaction UeP shows a cross section much larger than the reaction z/gC, so 
the number of events expected in future detectors and detected in SN1987 from the inelastic 
scattering is bigger than the elastic scattering. The direction of motion of the charged lepton 
on the incident neutrino is well preserved in i>e collisions. For reactions t^eP) the angular 
distribution is practically isotropic (there is a backward asymmetry). For elastic scattering 
ue it is expected that the distribution of the signal to be forward-peaked. Motivated by the 
anomaly of the directionality of detected events in Kamiokande-II and 1MB for SN1987A, it 
was proposed that some of the events detected were due to electronic neutrinos. However, 
such events had an unexpectedly high energy, but this energy should be smaller since the 
neutrino would carry out a part of it after the elastic scattering. Moreover, this process is 
extremely forward-peaked, especially for neutrinos with high energy, and hence such event 
candidates do not fit well the data from SN1987A (especially for 1MB) - for a recent analysis 
of events detected from SN1987A see [12]. 

For the next galactic supernova we expect a much larger amount of events coming from 
elastic scattering, then we consider that a good understanding of what happens to electronic 
neutrinos inside the supernova, in the context of neutrino oscillations, is crucial to understand 
the signal that will be detected. Needless to say that, for anti- neutrinos, the sample data 
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will be larger and with higher statistics the probability to find new phenomena will increase. 
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